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Abstract. We develop a general technique for proving convergence of re- 
peated quantum interactions to the solution of a quantum stochastic differen- 
tial equation. The wide applicability of the method is illustrated in a variety 
of examples. Our main theorem, which is based on the Trotter-Kato theorem, 
is not restricted to a specific noise model and does not require boundedness of 
the limit coefficients. 



1. Introduction 

It has been well established that the quantum stochastic equations introduced by 
Hudson and Parthasarathy [18] provide an essential tool in the theoretical descrip- 
tion of physical systems, especially those arising in quantum optics. The time 
evolution in these models is given by a unitary cocycle that solves a Hudson- 
Parthasarathy quantum stochastic differential equation (QSDE). These unitaries 
define a flow, which is a quantum Markov process in the sense of pQ, that rep- 
resents the Heisenberg time evolution of the observables of the physical system. 
Several authors have studied how quantum stochastic models can be obtained as 
a limit of fundamental models in quantum held theory [21 I15[ [9] . This provides a 
sound justification for using quantum stochastic models to describe physical sys- 
tems arising, e.g., in quantum optics. 

In contrast to the limit theorems for field theoretic models, which are non- 
Markovian and have continuous time parameter, it is natural to ask whether QSDEs 
can be obtained as a limit of discrete time quantum Markov chains. Classical coun- 
terparts of such results are ubiquitous in probability theory and there is a variety 
of motivations (to be discussed further below) to study such limits. The first re- 
sults on this topic date back to the work of Lindsay and Parthasarathy [26l [21] . 
In [26] it is shown that a particular class of repeated interaction models, where a 
physical system is coupled to a spin chain, converge in a very weak sense (in matrix 
elements) to the solution of a QSDE. A significant step forward was taken in [21] 
where the authors embed a chain of finite dimensional noise systems in the algebra 
of bounded operators on the Fock space and show strong convergence of the discrete 
flow to the flow obtained from a QSDE. Much later Attal and Pautrat [3] obtained 
similar results (in the special case of a spin chain) by showing that the discrete 
unitaries (rather than the flows) converge strongly to the solution of a QSDE. 

Independently from the previous work on the convergence of discrete chains, 
Holevo [161 117j has studied a very similar problem in his work on time-ordered 
exponentials in quantum stochastic calculus. The essence of Holcvo's approach is 
to define time-ordered stochastic exponentials as the limit of discrete interaction 
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models, where the role of the discrete noise is played by the increments of the field 
operators (the discrete noise is thus infinite dimensional in this setting, in contrast 
to the finite dimensional models considered by Lindsay and Parthasarathy). Despite 
the rather different motivation, the results of Holevo are strikingly similar to those 
obtained in the study of limits of discrete interaction models, as has been pointed 
out by Gough [14]. 

None of these results, however, are capable of dealing with the physically impor- 
tant case of limit QSDEs with unbounded initial coefficients (a typical setting in, 
e.g., quantum optics). This restriction is inherent to the techniques used to prove 
these results, which rely on Dyson series expansions and require boundedness of the 
coefficients. Moreover, each of the above results has been proved separately in its 
own setting, while the similarity between these results strongly suggests that they 
should be unified within a common framework. 

The purpose of this paper is to introduce a general technique for proving conver- 
gence of a sequence of discrete quantum Markov chains to the solution of a QSDE. 
Our approach does not rely on a Dyson series expansion, but instead employs a 
form of the Trotter-Kato theorem. This allows us to deal with unbounded coeffi- 
cients in a natural and transparent manner. In the simpler case where the limit 
coefficients and/or discrete noises are bounded, we obtain many of the previous 
results as special cases of our general theorems. Moreover, the specific functional 
form of the limiting coefficients, obtained in [3l ll61fT4] by identifying a power series 
obtained from the Dyson expansion, is effectively demystified: we will see that it it 
is an immediate consequence of the scaling of the noise operators. 

Our motivation for this work is twofold. First, the convergence of discrete quan- 
tum Markov chains to continuous ones is a fundamental problem in quantum prob- 
ability theory. In classical probability theory such problems have been investigated 
for many decades, and the theory has culminated in the well known work of Stroock 
and Varadhan [30] section 11.2], Ethier and Kurtz [10] and Kushner [20], among 
others. A similar systematic investigation was hitherto lacking the the quantum 
probability literature. This paper presents one attempt to unify and extend the 
existing results in this direction. 

Second, we are motivated by practical problems in which one is specifically in- 
terested in the convergence of discrete to continuous models. For example, certain 
laboratory experiments, e.g., atomic beam experiments with a large flux of atoms, 
can be approximately modelled by quantum stochastic equations in the limit of 
large flux. Another application of independent interest is the development of nu- 
merical methods for quantum stochastic models. To perform tractable numerical 
simulations one is often forced to discretize, particularly in dynamical optimiza- 
tion problems which appear in the emerging field of quantum engineering [6] , and 
convergence of the discretized approximations is a challenging topic. We were mo- 
tivated in particular by the problem of discretizing linear quantum systems, which 
play a special role in linear systems theory [T9l [25] [24], but are not covered by 
previous results as both the noise and the initial system are necessarily unbounded. 

As compared to previous results, our approach is closest to the original method 
of Lindsay and Parathasarathy [21] , The simple uniform convergence result [HI 
proposition 3.3] is replaced in our setting by a variant of the Trotter-Kato theorem 
[10] . which allows us to deal with the analytic complications inherent to the case 
of unbounded coefficients. We also work directly with the unitary evolution rather 
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than with the flow. The Trotter-Kato theorem allows us to obtain convergence by 
studying generators, and exploits in a fundamental way the Markov property of 
both the approximate and the limit evolutions. 

Techniques to obtain convergence for QSDEs by studying generators were intro- 
duced by Fagnola [12] and Chebotarev [8] (using resolvents) and by Lindsay and 
Wills [221 [23] (using the Trotter-Kato theorem). We have previously applied related 
techniques to obtain general results on singular perturbation problems for quantum 
stochastic models [5[ [7] . The application of techniques of this type to obtain the 
convergence of discrete quantum models is new. 

The remainder of this paper is organized as follows. In section [2] we introduce 
the class of discrete interaction models and limit models which will be of interest 
throughout the paper, and we state our main results. The main theorem is a 
generalization of the Trotter-Kato theorem to quantum stochastic models, and is 
generally applicable. We also introduce a more restricted family of discrete models 
for which the conditions of this result can be verified explicitly. Section [3] develops 
a variety of known and new examples using our results. Finally, section[4]is devoted 
to the proofs of our results. 

2. Main Results 

In the following subsections, we first define the class of models that we will 
consider and introduce the necessary assumptions. This is followed by the statement 
of our main results. The proofs of our main results are contained in section [4] below. 

2.1. The limit model. Throughout this paper we let TL, the initial space, be 
a separable (complex) Hilbert space. We denote by T — T s (L 2 (R + ;C n )) the 
symmetric Fock space with multiplicity n £ N (i.e., the one-particle space is 
C" <g> L 2 (R+) = L 2 (M+;C n )), and by e(/), / G i 2 (R+;C n ) the exponential vec- 
tors in T . The annhiliation, creation and gauge processes on J 7 , as well as their 
ampliations to TL ® T, will be denoted as A\ , and A 1 / , respectively (the chan- 
nel indices are relative to the canonical basis of C"). Moreover, we will fix once 
and for all a dense domain T> C TL and a dense domain of exponential vectors 
£ = span{e(/) : / G &} C J 7 , where 6 C L 2 (K+;C n ) (~l L 1 ~(K+;C Tl ) is an ad- 
missible subspace in the sense of Hudson-Parthasarathy [18] which is presumed to 
contain at least all simple functions. An introduction to these concepts (using a 
similar notation to the one used here) can be found in [4] , and we refer to [HI [22] 
for a detailed description of quantum stochastic calculus. 

Consider a quantum stochastic differential equation of the form 

{n 7i n 

J2 (Nij - Sij) dA? + M i dA t + ^L l dAl+K dt 
i,j=l i—1 i—1 

where Uq = I and the quantum stochastic integrals are defined relative to the 
domain T> ® £ (for simplicity, we will use the same notation for operators on TL 
or on T and for their ampliations to TL® T). Under certain conditions to be 
introduced below, the solution of this equation describes the time evolution of 
quantum stochastic models such as those used in quantum optics. The purpose of 
this paper is to prove that the solution of this equation may be approximated by 
appropriately chosen discrete interaction models. The approximating models will 
be defined in the following subsection. 
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Denote by 6 t : L 2 ([t, oo[; C") -> i 2 (R+;C") the canonical shift 6 t f(s) = f(t+s), 
and by Ot : Tu — > T its second quantization (here T = T t \ ® denotes the 
usual continuous tensor product decomposition). Recall that an adapted process 
{Ut : t > 0} on TL ® T is called a unitary cocycle if Ut is unitary for all t > 0, 
i h- > ZTj is strongly continuous and U s +t = U S (I ® 0* s Ut<d s ), where 7 ® 9*[/t0 s is 
viewed as an operator on J 7 ^ ® (7i .Fr s ) = TL ® T. The following condition will 
always be presumed to be in force. 

Condition 1. TTie operators K, Li, Mi, and Nij, defined on the domain T>, are 
such that the Hudson-Parthasarathy equation (Qp possesses a unique solution {Ut : 
t > 0} which extends to a unitary cocycle onTL® T '. 

Remark 1. When the coefficients K, Li, Mi, iVy are bounded, it is well known 
[18] that condition [T] holds true if and only if the following algebraic relations are 
satisfied: 

n n 

K + K* = -J2 LiL*, Mi = -J2 NijL*, 

n n 
3=1 3=1 

Remarks on the verification of condition [T] in the unbounded case can be found in 
0. 

Remark 2. We have chosen the left Hudson-Parthasarathy equation fl}, rather 
than the more familiar right equation where the solution is placed to the right 
of the coefficients. This means that the Schrodinger evolution of a state vector 
i/j € TL (8 T is given by U^tp, etc. The reason for this choice is that for equations 
with unbounded coefficients, it is generally much easier to prove the existence of a 
unique cocycle solution for the left equation than for the right equation (see, e.g., 
[TTJ[52]). As we will ultimately prove convergence of discrete evolutions to t/ t *, there 
is no loss of generality in working with the more tractable left equations. If we wish 
to begin with a well defined right equation (as is more natural when the coefficients 
are bounded) , our results can be immediately applied to the Hudson-Parthasarathy 
equation for its adjoint. 

2.2. The discrete approximations. A discrete interaction model describes the 
repeated interaction of an initial system with independent copies of an external 
noise source. Given an initial Hilbert space TL' and a noise Hilbert space K! , a 
single interaction is described by a unitary operator R! on TL' ® K! '. To describe 
repeated interactions, we work on the Hilbert space TL' <E) ® N K,' on which we define 
the natural isomorphism 

S fc : (/cT (fe - 1} ® H' ® (g) K' -»■ W ® (g) K! 

N N 

as Ek{ipk-i] <8 £ <8> tp[k) = £, <8 (ipk-i] ® ^[k)- We now define recursively R' Q = I and 

Kit ® rpk-i] ® ipk ® ip[k+i) = K-i^k(ip k -i] ® R'(£ O ipk) ® -0[fc+i) 

for k e N, where ip k _ 1} ® Vfc ® ^[fc+i € (/C')® (fc_1) <8> JC' <8> ® N /C' ® N JC' and 
£ € W. This is precisely the discrete counterpart of a unitary cocycle. Note that 
the order of multiplication in the recursion for R' k matches the choice of the left 



DISCRETE APPROXIMATION OF QUANTUM STOCHASTIC MODELS 5 

Hudson-Parthasarathy equation above (i.e., R'^tp is the Schrodinger evolution, see 
remark [2J . 

The purpose of this paper is to prove that the solution of the Hudson-Parthasarathy 
equation fl} can be approximated by a sequence of discrete interaction models with 
decreasing time step. In order to study this problem, we will embed our discrete 
interaction models in the limit Hilbert space Ti®J-. This allows us to prove strong 
convergence of the embedded discrete cocycles to the solution of equation ^ . The 
precise way in which the embedding is done does not affect the proof of our main 
result; we therefore proceed in a general fashion by defining a fixed but otherwise ar- 
bitrary sequence of discrete interaction models which are already embedded into the 
limit Hilbert space (this avoids, without any loss of generality, the notational bur- 
den of introducing separate Hilbert spaces and embedding maps for every discrete 
approximation). As a special case of this general construction, we will introduce in 
section below an interesting class of discrete interaction models for which the 
embedding is made explicit. 

We proceed to introduce the embedded discrete interaction models. For every k £ 
N, we will define a discrete interaction model with time stepQ 2~ k (the ultimate goal 
being to take the limit as k — > oo). By its continuous tensor product property, the 
Fock space is isomorphic to T = (^) N T 2 -k-\ , where each component T 2 -k\ represents 
a consecutive time slice of length 2~ k . Our discrete interaction models will be 
embedded as repeated interactions with consecutive time slices of the field. Let us 
make this precise: the following notations will be used throughout the paper. Let 
Ti k C Ti be the initial Hilbert space and let K. k C J- 2 -k] be the noise Hilbert space 
for the interaction model with time step 2~ k . We will write 

N N 

We now introduce a unitary operator R k on Ti k ®JC k which describes the interaction 
in a single time step, and we extend R k to Ti ® T by setting R k ip = ip for ip £ 
(Ti k ®IC k (g) J" [2 -0 X - We now define recursively R k = I for t £ [0, 2~ fc [ and 

R k = R k (g _ 1)2 - k {I ®<d* (l _ 1)2 - k R k <d { ^ 1)2 - k ) forte [£2" fc , + l)2" fc [, feN, 

where I ®Q* t2 _ k R k @i 2 -h is viewed as an operator on T t2 -k^®(Ti®Ty l2 -k) = Ti®F. 
In words, the discrete evolution R k is an adapted unitary process which is piecewise 
constant on consecutive intervals of length 2~ k , and an interaction between the 
initial system and the next field slice occurs at the beginning of every interval. Our 
goal is to prove that R k converges as k — > oo to the solution Ut of equation |T]) . 

Remark 3. This model coincides precisely with the discrete interaction model 
defined above if we choose H' = H k , K' = K k , R' = R k , and R' t = R k 2 _ n . 

2.3. A general limit theorem. Before we proceed to the statement of our main 
result, we introduce certain families of semigroups which will play a central role in 
our approach. 

Lemma 1. For k £ N and *p,<p £ T 2 -k\, define R k '^^ • Ti — * Ti such that 
(u : R k ^v) = HVll^lMrV® i>,R k v®<p) Vu,v£ Ti. 



This choice is only made to keep our notation manageable, and is by no means a restriction 
of the method. 
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Then R^v is a contraction on 7i and 

(u, {R k ^) L t2 %) = U\\- N2 " y\\- N2 " (u ® i>® N2k ,R k t v® v® N2 "} 
for all u,v EH and t G [0,N], N G N. 

Proof. This follows directly from the unitarity of R k and the definition of R k . □ 

The following counterpart for the limit equation ([T]) is proved in [7J lemma 1] . 
Lemma 2. For a, (3 G C n , define T t Q/3 :TL^TL such that 



(u, T t al3 v) = e -tt a \ 2+ W 2 W 2 (u ® e(aZ [0 ,t]), U t v® e{f3I m )) Vu,veH, t> 0. 

Then T"' 3 is a strongly continuous contraction semigroup on Ti., and the generate 
J£ a P of this semigroup satisfies Dom(J>f a ' 3 ) D V such that for u G T> 



( n n n I |2 _|_ I o|2 

Y aSNijPj + J2 a * M i + E L ^ + K ~ 2 
i.j — 1 i—1 i—1 

The reason to focus on the semigroups associated to our models is that we 
will seek conditions for convergence of the discrete approximations in terms of the 
generators j£f Q/3 . As the latter are expressed directly in terms of the coefficients of 
the limit equation JT]), such conditions can typically be verified in a straightforward 
manner and do not require us to work directly with the solution Ut of that equation. 

The following is the main result of this paper. The theorem bears strong re- 
semblance to the Trotter-Kato theorem for contraction semigroups, and the latter 
does indeed form the foundation of the proof. The proof of the theorem is given in 
section [4] 

Theorem 1. Assume that conditional holds, and let T> a/3 C Dom(Jzf a ' 3 ) be a core 
for ££ a P , a, (3 G C". Then the following conditions are equivalent. 

(1) For all a,/3 G C™, u G V afj there exist u k G H and ip k , tp k G F 2 -k] so that 

u (-0) ► e(aJ[ ,i]), VP) > e(/?/[o,i]), 

and 

(2) .For every T < oo and ip £ H. ® J 7 

lim sup ||i? t fe * V — = 0. 

fc^oo o<t<T 

Remark 4. As was pointed out to us by a referee, this theorem could also be 
stated outside the framework of quantum stochastic differential equations. Indeed, 
a careful reading of the proof reveals that one may dispose of condition [1] entirely 
and assume only that Ut is any unitary cocycle on H. (g> T . The theorem is the 
most powerful, however, when the generators «5f a ^ admit an explicit expression in 
terms of the parameters of the limit model. This is the case, in particular, when 
condition [1] is satisfied and T> is a core for all Jz?"^, a,f3 G C. We can then choose 
D a P C T>, with the important consequence that this puts the explicit expression 
for j£? a in lemma [5] at our disposal. This will be the case in all our examples. 
Typically existence and uniqueness proofs for the solution of |T]) already imply that 
I? is a core for ££ a $ ; see, e.g., [11] [22] and [3 remark 4] for further comments. 
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Remark 5. The assumption in condition [T] that Ut is a unitary cocycle can be 
weakened somewhat. In the absence of this assumption, one may still prove the 
implication 1=>2 provided that strong convergence uniformly on compact intervals 
is replaced by weak convergence of R k to Ut for every time t. We have chosen to 
concentrate on the unitary case, as it is the relevant one for physical applications 
and admits a much stronger result. 

2.4. A class of discrete interactions. Theorem[T]is a very general result which 
allows us to infer the convergence of a sequence of discrete interaction models 
to the solution of the Hudson-Parthasarathy equation [TJ The verification of the 
conditions of the theorem requires additional work, however, and the form of the 
limit coefficients depends on the choice of the discrete interactions R k . In this 
section we introduce a special class of discrete interaction models (with Ti k = Ti) 
for which the conditions of theorem [T] can be verified explicitly. In particular, we 
obtain explicit expressions for the limit coefficients. It should be noted that this 
class of discrete interaction models is physically natural, and we will encounter 
several concrete examples in section |3l 

Let IC be a fixed Hilbert space and suppose that we are given a sequence of 
bounded operators n k : IC — ► T 2 -^ which are partially isometric in the following 
sense: 

■K k *K k = I Kj n k ir k * — P K k for all k, where )C k := ran ir k . 

Here P K k is the orthogonal projection onto IC k and Ik. is the identity on IC. The same 
space K. will play the role of the noise Hilbert space for every discrete interaction 
model R k (fc £ N) after being isometrically embedded into the limit Hilbert space 
by the embedding maps ir k . As we will see, the choice to work with a fixed noise 
space prior to embedding is convenient due to the fact that the limit coefficients 
can be expressed in terms of the matrix elements of certain operators on Ti ® IC. 

We now define the discrete models. Let F±, . . . , Ft, G\, . . . , G m , and H±, . . . , H r 
(£, m,r£N) be bounded self-adjoint operators on Ti, and let Ai, . . . , A^, fii, . . . , /i m 
and ui, . . . , v T be (not necessarily bounded) self-adjoint operators on /C. Define 

i m r 

H k = 2 k Y J Fj ® <\j + 2 fe / 2 G i ®H+Y1 H i v i 

2=1 2=1 2=1 

on H®(f)j Dom(Aj) n p|j Dom( ft ) n fij Dom(i/ i )) := V . 

Condition 2. The following are presumed to hold: 

(1) H is essentially self-adjoint for every k. 

(2) F := F\ ® \\ + • • • + Fg (8 \i is essentially self-adjoint on the domain T>q. 

(3) There is a family of orthonormal vectors xoi ■ ■ ■ j Xn G IC such that for all 
a e C n , the vector x k (a) ■= xo+2^ k ^ 2 YTj=i a jXj satisfies (ir k x k (a))® 2k — » 
e(a/[o,i])- 

(4) xo € Dom(yj) for all j. 

(5) xo G Dom(ft) and (xa,VjXo) = for all j. 

(6) Xo € Dom(Aj) and XjXo = for all j. 

We subsequently identify H k and F with their closures. We now define the 
discrete interaction unitary R k :— Tr k e lH 2 n k * on H ® IC k , and extend to Ti (£> T 
as usual. 
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To state the convergence result for this class of discrete interaction models we 
must introduce the corresponding limit coefficients Nj, Mi, Li and K in equation 
([1}. This is what we turn to presently. Define the bounded continuous functions 
/, g : R -> C as 

f(x) = , g(x) = - 2 . 

We must also define the bounded operators Wij , Xf, Yf , Z pq (i,j — 1, . . . , m, p, q = 
1, . . . , n) on TL as follows: for every u, v £ H, 

(u,WijV) = (u<E> ^xo,g(F)v<g) HjXo), 

(u,Xfv) = {u®x P J(F)v®fiiXo), 

(u,Y?v) = (u<g> /J-iXo,f(F)v®x P ), 

(u,Z pq v) = {u® X p,e lF v®x q )- 

It is evident that these operators are well defined provided that condition [2] is 
assumed to hold. We now define the limit coefficients as follows. 

rn m 

N pq = Z™, M p = X?Gi, L P =J2 GiYf, 
i=l i=l 
r m 

K = iJ2Hj (xo, v jX o) + G i w H G r 

3=1 »i3=l 

The following is the main result of this section. 

Theorem 2. Suppose that condition^ holds and that condition^ holds for the limit 
coefficients N pq , M p , L p , K as defined above (with T> — TL). Then 

lim sup \\B$*ip - U?ip\\ = for all ip £ TL ® T , T < oo. 

fe^oo o<t<T 

Remark 6. We have not sought to develop this result under the most general 
conditions possible. In particular, the boundedness of Fj,Gj,Hj can certainly be 
relaxed if appropriate domain assumptions are introduced, and the proof of the 
present result is then readily extended. We have chosen to restrict to the bounded 
case as the treatment of this case is particularly transparent, and we have not found 
one single choice of domain conditions in the unbounded case which covers all ex- 
amples of interest in that setting (particularly straightforward extensions can be 
found when either Fj = for all j, or when only the Gj are bounded). An illus- 
trative example with unbounded coefficients is developed in section [3] by appealing 
directly to theorem [TJ 

One might worry about the well-posedness of theorem [5] as the limit coefficients 
Npq, M p , L p , K depend on our choice for %q, ■ ■ ■ > Xn- For completeness, we provide 
the following simple lemma whose proof can be found in section 2] 

Lemma 3. Define x k ( a ) as ^ n condition^ (which we presume to be in force). 
Suppose that XOt ■ ■ ■ ->Xn € is another orthonormal family such that for every 
a e C" 

n 

X k (a) := xo + 2" fc / 2 £ a jXj satisfies (ir k X k (a))® 2 " e(aI [0A] ). 

3=1 

Then there is a <f> £ M such that Xj = e Xj f or a ^ j ■ 
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Finally, it should be noted that condition [T] imposes stronger assumptions on 
the noise operators fj,j and Xj than is evident from condition O The following 
lemma provides explicit assumptions which guarantee that condition [1] holds for 
the coefficients N pq , M p , L p , K as defined in this section. The proof is given in 
section |4j 

Lemma 4. Define S = span{%i, . . . ,Xn}, and suppose that fJ-jXo G <5 for all j. 
Suppose moreover that S C Dom(Aj) and that XjS C S for all j. Then the limit 
coefficients N pq , M p , L p , K satisfy condition^ 

Remark 7. It is not difficult to construct examples of discrete interaction models 
where the limit coefficients N pq , M p , L p , K do not satisfy the Hudson-Parthasarathy 
conditions in remark Q] In this case, our proofs may be modified to show that R k 
nonetheless converges weakly to the solution U t of equation (JTJ) with coefficients 
N pq , M p , L p , K as defined in this section. However, in this case the limit evolution 
Ut will not be unitary, so that the physical relevance of such a result is rather 
limited. 

3. Examples 

We illustrate our results using four examples. The first two examples reproduce 
results of Attal and Pautrat [3 and Holevo [16j . The third example, that of a 
linear quantum system, possesses both unbounded noise operators and unbounded 
initial coefficients. Finally, the fourth example shows that we may simultaneously 
approximate the noise and the initial coefficients, as is often useful in numerical 
simulations. 

Until further notice TL is a fixed Hilbert space. Note that for sake of example, 
all our models live in the Fock space with unit multiplicity n = 1. The extension 
to multiple channels is entirely straightforward and leads only to complication of 
the notation. 

3.1. Spin chain models. Define JC = C 2 and denote the canonical basis of C 2 as 
(XOjXi)- The noise Hilbert space is thus that of a single spin. We embed the spin 
into the Fock space by defining the embedding map ir k : JC — > ^2- fe ] through 

* k Xo = e(0) = 10 00, ir k X i = 2 k ' 2 I [Q>2 - k] ® 

p=l p=2 

(here /[o,2- fc ] is the indicator function on the interval [0, 2~ fe ]). We now define 
bounded self-adjoint operators A, [i\ and /i2 on IC as matrices with respect to the 
canonical basis: 

A=(J [J), /ii=^=(j J), M2 = ^=(° 7). 

Let F, Gi, G2 and H be arbitrary bounded self-adjoint operators on Tt, and let 
Hk be a self-adjoint operator on IC. Clearly the operator H k defined on Ti. <S> IC by 

H k = 2 fe F <g> A + 2 k ' 2 (Gi ® m + G 2 ® /i 2 ) + H ® I + I ® H K 

is self-adjoint for every k, and the domain assumptions in condition [2] are trivially 
satisfied. Furthermore, we have Axo = and (xa^iXo) = 0- 
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For every a € C, we have defined x k ( a ) — Xo + 2~ fe / 2 axi- Note that 

\\e(al [0ll] ) - (n k X k (<x))® 2k \\ 2 = \\e(al m - k] )f +1 + \\n k X k (a)\f +1 
- «uI [0t2 - k] ),iT k x k (a)) 2k - {TT k X k {a),e{aI [0t2 - k] )f k 

= el Q l 2 - (1 + \a\ 2 2- k f 0. 

We conclude that (Tr k X k 

( ))»2* ±^ e ( a J [0jl] ). We have now verified all parts of 
condition^ Moreover, the coefficients Nu, Mi, Li and K are easily computed: 

%F r %F j 

N n = e iF , M 1 = ^-—(Gi-iG 2 ), U = (Gi + iG 2 ) 

r r 

K = iH + i(xo,H K XQ) I + (Gi + *Ga) ^ ~ *f ~ 1 (Gi - »G 2 ). 

Note that these coefficients satisfy condition Q] by virtue of remark [TJ It follows 
from theorem [2] that the repeated interaction R k * obtained from the Hamiltonian 
H k converges strongly to the solution C7 t * of equation ([1]) with these coefficients, 
uniformly on compact time intervals. This result corresponds to [3J theorem 19]. 

3.2. Time-ordered exponentials. Let F, G, H be bounded operators on Tt with 
F, H self-adjoint, and define the essentially sclf-adjoint operators 

AM k := F (A te - A t< _J + G (A\ t - A\ t _J + G* (A u - A t ,_J + 2~ fc 

on Tt®£, where ti — £2~ k . Holevo [16] defines time-ordered stochastic exponentials 
as the following strong limits on Ti. ® T: 



exp 



i f (FdA s + GdAl + G*dA s + H 
Jo 



ds) 



:= s-hm e L«*J • • • e lAAil . 

k — >oo 



Evidently this definition can be interpreted as the limit of a sequence of discrete 
interaction models, and the limit process does indeed solve an equation of the form 
([1]). In this example we develop this idea by applying theorem [2] 

Define K. = T-n = r s (L 2 ([0, 1])), and choose the orthonormal vectors xo and Xi 

as 

oo oo 

Xo = e(o) = i®0o, » = oe/ M e0o. 

We define embedding maps ir k : K. — » T 2 -^\ by specifying their action on the (total) 
family of exponential vectors e(/) € T\\, f € L 2 ([0, 1]) as follows: 

T fe e(/) = e(/ fe ), f*{x) = 2 k / 2 f(2 k x) 

(note that f k G L 2 ([0, 2~ k ])). It is easily verified that ir k is a unitary map for every 
k. We define self-adjoint operators /Ui, and A on K, as 

^i=Ai + At, fj, 2 = i(Ai- A\), A = Ai, 

where j4i,Aj and Ai are the standard Hudson-Parthasarathy noises evaluated at 
time * = 1. Define Gi = (G + G*)/2 and G 2 = i(G - G*)/2, and set 

H k := 2 k n k *AM k T: k = 2 k F ® A + 2 k/2 (Gi <g> //i + G 2 ® {i 2 ) + H ® I . 

It is well known that and F = F ® A are essentially self-adjoint, and it is easily 
verified that (xo,MiXo) = and Axo = 0. Moreover, as 7r fe xo and 7r fe xi coincide 
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with their counterparts in the previous example, we have verified that condition [2] 
holds. The coefficients Nn, Mi, L\ and K are now easily computed: 

p iF — T p iF — T 

N n = e lF , M l = -— -- G, L x = G* — — , 

r r 

e lF - I - iF 

K = iH + G* -i— -- G - 

T z 

These coefficients satisfy condition [T] by virtue of remark[TJ It follows from theorem 
[5]that the time ordered exponential defined above coincides with the adjoint solution 
C/ t * of equation |T]) with these coefficients. This agrees with [THl corollary 1]. 

3.3. Linear quantum systems. Let TL = £ 2 (Z+) and denote by (rjk, k G Z+) the 
canonical basis in £ 2 (Z+). We also choose the domain V = span{j7& : k G Z+} C TL 
of finite particle vectors. On T> we define the operators 

arjk = Vkrik-i, a* rf k = \Jk + 1 r] k +i, q = a + a* 1 p = i(a-a*). 

Note that a is the annihilation operator and a* is the creation operator, while q 
and p are the position and momentum operators, respectively. 

A linear quantum system is a quantum stochastic differential equation of the 
form (p} whose coefficients take the following form on T>: 

N\i=I, Mi = mp + m' q, Li — — m* p — to'* q, K = iH + — L\M\, 
H = kip 2 + k 2 (pq + qp) + k 3 q 2 + k 4 p + k 5 q + k 6 I, 

where to, to' G C and ki,...,kQ £ K. Physically, a linear quantum system is a 
model whose Hamiltonian is quadratic in position and momentum and whose noise 
coefficients are linear in position and momentum. At least formally, one may easily 
verify that the Hudson-Parthasarathy conditions of remark[T]are satisfied, but the 
coefficients are unbounded in this case. That condition [T] is satisfied in this setting 
is proved in [11) . 

Linear quantum systems possess various special properties: for example, the 
adjoint solution ?7 t * of equation dXJ) leaves the family of Gaussian states in TL ® T 
invariant, the Hcisenberg evolution of the observables {q,p) has an explicit solution, 
and the quantum filtering problem for ([T]) has a finite-dimensional realization (the 
Kalman filter). Because of these and other properties, the linear quantum systems 
play a special role in quantum engineering as they admit particularly tractable 
methods for control synthesis and signal analysis [TH1 [551 US] • In these applications 
it could be of significant interest to work with discrete time approximations (e.g., for 
the purpose of digital signal processing) , but it is important to seek approximations 
which preserve the linear systems properties of these models. This is easily done, 
but we necessarily obtain discrete approximations where both the initial system 
coefficients and the discrete noises are unbounded. In this example we will prove 
the convergence of such unbounded discrete approximations by appealing directly 
to our main theorem [1] 

We begin, however, by setting up our discrete models as in section 12.41 Let 
K. = ^ 2 (Z + ), and choose an embedding ir k : K. — > T- 2 -^ by setting 

oo i— 1 oo 

Ao = 1 6 00, n k m = 00© (2 fe / 2 I [0 , 2 - k] ) m © (£ G N). 
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If we define r] k (a) := tjq + 2 k / 2 oerii for a £ C, then we find precisely as in the 
previous examples that (w k r\ k (a))® 2 * — ► e(aJ[o ; i]) as k — > oo. 
Let us now define on 2? ® 2? the operators 

H k = H®I -i (Mi a* + Li <g> a) 2 fe/2 . 

One may verify that H k is symmetric and that T> <g> T> is a domain of analytic vectors 
for iJ fe , so that in particular H k is essentially self-adjoint for every k [2"5l section 
X.6]. We will subsequently identify these operators with their closures. Note that 
the Hamiltonian H k is quadratic in the family of position and momentum operators 
of the initial system and of the discrete noise; therefore the discrete interaction 
model is itself a (discrete time) linear system, and it therefore possesses all the 
associated desirable properties. 

We now define the discrete interaction unitary R k on TL^T from the Hamiltonian 
H k as in section HOI We will use theorem [T] to prove that 

lim sup \\R k *ip - U^W = for all ip £ H ® T, T < oo, 

fc-»oo o<t<T 

where Ut is the solution of equation (JT|) with the coefficients Nu, Mi, L\, K defined 
above. Note that theorem[5]does not apply as the initial coefficients are unbounded, 
but we may essentially repeat the proof of that theorem with minor modifications 
to obtain the present result. To this end, we begin by noting that I? is a core for 
Jzf"' 3 (a, (3 £ C) by the analytic vector theorem (see remark 4]). Let us fix 
a, f3 £ C, and define ip k = Tr k rj k (a) and ip k — ir k r] k (j3). By theorem[IJ it suffices to 
prove that 

|| 2 fc (jR fc;VV _ j) u _ j^a/Juji o 

for every u £ T>. We now proceed as follows. Fix u £ T>, and note that using the 
trivial identities e lx = 1 + f(x) x = 1 + ix + g(x) x 2 we can write 

e lH " 2 ~ k u(g> m = {I + 2- k f(H k 2- k ) H k ) u(g>g e 

= (I + iH k 2- k + 2- 2k g{H k 2- k ) {H k ) 2 ) u <g> m . 

Here we have used the spectral theorem and the fact that u <g> r/i £ Dom((H k ) p ) for 
every £,p (see the proof of theorem [2] for a more precise argument). Therefore 

(v <g> r 1 k (a),2 k (e lHk2 ~ k - I) u ® rj k ((3)) = i(v, Hu) 

- (v <g> r]o, g(H k 2~ k ) (Mi ® a* + L% ® a) Mm <g> 771) 

- i a* (v <g> 77!, f(H k 2- k ) M lU <g> 771) 

- i (u <8> i]o,f(H k 2~ k ) (Mi (g) a* + ii (g) a) u (g) 771) /3 
+ a* (7; <g> ft, ( e ^ fc2 ~ fc - /)« ® rri) /? + 0(||«|| 2- fe / 2 ). 

A straightforward computation shows that 

sup \{v®r) k {a),2 k {e iH " 2 ~ k - I)u®rf{p)) - (v,(a* M x + Li{3 + K)u)\ 0, 

vt=n 
\\v\\<l 

where we have used that f(H k 2- k ) -> iJ, g{H k 2- k ) -> -(1/2)7 and e^* 2 "* _> J 
strongly as fc — > 00 by theorems VIII. 20 and VIII. 25]. As in the proof of 
theorem [2] below, it follows readily that \\2 k (R k ^ k - I)u - 5? af) u\\ -> 0. 
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3.4. Finite dimensional approximations. In the previous examples we have 
approximated the quantum stochastic differential equation ([1]) by constructing dis- 
crete interaction models whose interaction unitary R k lives on the Hilbert space 
TL ® K. k . Even though the Fock space T is infinite dimensional, we have seen 
that we may choose finite-dimensional discrete noise spaces as simple as JC k = C 2 . 
In numerical applications, however, we typically wish to go one step further and 
approximate also the initial space TL (which is often infinite dimensional) by a fi- 
nite dimensional space TL k . The discrete interaction models then live entirely on 
the finite dimensional Hilbert spaces TL k <g> K. k , as is desirable for numerical imple- 
mentation. We would like to establish that these discrete models converge to the 
solution of the limit equation |T]) when we simultaneously let the time step go to 
zero and let the initial space dimension go to infinity. We will now show in a toy 
example that this problem fits into the setting of theorem [TJ 

We will discretize the noise essentially as in the example of section |3"7T1 but let 
us directly embed the discrete models (rather than work with the embedding maps 
TT k ) to simplify the notation. For every k £ N, define the vectors XoiXi € 3~2- k \ as 

oo oo 

Xo = l©0O- Xi =O®2 fe / 2 / [Oi2 - fc] ©0O. 

p=l p=2 

Moreover, let TL be an infinite dimensional separable Hilbert space, and fix an 
orthonormal basis {r/k ■ k € N} C TL. For the ktYi discrete interaction model we will 
choose the initial Hilbert space Ti k := span{?7i, . . . ,i] k } C TL and the noise Hilbert 
space IC k := span{xo,Xi} C T^-hy Thus TL k JC k is indeed finite dimensional. 

Let H and M be bounded operators on TL where H is self-adjoint. In this toy 
example, we will be interested in approximating the solution of the limit equation 

dU t = U t {M dA\ - M* dA t +iHdt- \M*M dt}, 

which satisfies condition [1] by virtue of remark [T] We claim that this can be done 
by choosing the discrete interaction unitaries 

R k = exp{2-' £ / 2 (P k MP k ®b* k - P k M*P k ®b k )+i 2~ k P k HP k <g> I}, 

where P k denotes the orhogonal projection onto TL k and the noise operator b k is 
defined by setting b k Xo = 0, b k Xi = X§' an d b k ip = for ip _L K. k . 

We simply verify the conditions of theorem[TJ As all the coefficients are bounded, 
we may choose V a P — T> — TL. Fix a, (3 € C and u £ TL, and let us define 
^ fe = Xo + 2 _,£ / 2 axi an d f k = Xq + 2~ fe / 2 /3xi As in the previous examples we 
have 

WT 2 " — e{oJ m ), &T 2k — e(fll m ). 
As P k u — > u as k — > oo, it suffices to verify that 

2 k {R k^^ _ I)PkU ±^ ^ a p u _ 

What remains is again essentially the same computation as in the previous example 
and as in the proof of theorem [2j We leave the details to the reader. 

4. Proofs 



4.1. Proof of theorem [U The proof of theorem Q] is based on a version of the 
Trotter-Kato theorem for contraction semigroups due to Kurtz. We cite here from 
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[TO], theorem 1.6.5] a special case of this result in the form which will be convenient 
in the following. 

Theorem 3 (Kurtz). Let TL be a fixed Hilbert space. For k £ N, let T k be a linear 
contraction on TL and let T t be a strongly continuous contraction semigroup on TL 
with generator Jz? . Let T> be a core for . Then the following are equivalent. 

(1) For every u € T>, there exist u k £ TL such that 

ufej fc- ! oo^ 2 k (T k -I)u k ±^>_S? U . 

(2) For every ip £ TL and t < oo 

Km ||(T fc )L t2fe Jv--T^|| =0. 

k — >oo 

(3) For every ip gH and t < oo 

lim sup||(T fc )L s2fc J^-7;V>|| =0. 

Armed with this result, we may now proceed to prove theorem [1] We will first 
prove the forward direction, and we subsequently consider the converse implication. 

Theorem^ 1=>2. Let us restrict our attention to the interval [0, N] with N £ N. 
It suffices to prove that convergence holds uniformly on [0, N] for any N. We may 
therefore restrict the Hilbert space to TL (g> Tm\ , which we do from now on. 

First, let a, (3 £ C n and ip k , ip k £ 3~2~ k ] as m the statement of the theorem. Then 

(u® {il> k )® N2 \R k t v® (<p k )® N2k ) = U k \\ N2k \\ V k \\ N2k (u,(R k ^ k ^)L t2k iv) 

— e (|Q|2+l/3|2)A,/2 <", Tfv) ={u® e(aI [0 , N] ), U t v ® e (/?/ [0>J v])> 

for any u,v £ TL and t G [0, iV] by lemma Q] and theorem [3l Similarly, we can 
establish the following. Let = to < t\ < ■ ■ ■ < t m < t m +l — N (m 6 N) be a 
dyadic rational partition of [0, N], i.e., tj — £j2~ k for some k' £ N with ij £ N for 
all j = 1, ... , m. Let ao, . . . , a m , . . . , (3 m £ C™, and choose tp k , (p k £ T 2 -^ such 
that 

i^r 2 " — e (Mo,n), ^r 2 " — e(/3,i [0>1] ). 

Let /,g 6 L 2 ([0, iV]; C") be simple functions with f(s) — ctj and g(s) = (3j for 
s £ [tj, tj+\[, and define for all k > k' the vectors 

^ = (^ k )^ 2k - k ' ® ( ^)8(fe-«Oa*-*' 8 . .. ® ( ^,fc ) ®(wa*'-/ m )a»-* , j 

Then — ► e(/) and — > e(g) as fc — > oo, and it is not difficult to verify (using 
the cocycle property of Ut) that for t £ [tj,tj+i[ 

{u®e{f),U t v®e{g)) = ||e(/)|| \\e(g)\\ (u,T t a f°T^ t \ ■ ■ -T^v). 

We therefore obtain for u,v £ TL and t £ [tj,tj+i[ 

(u®Vk, = ||^fc||||^fc||<« J (i2 fc ^) tia ^--(fl^^) L(t - t ^ 2 * J t;) 

^ || e(/) || || e(ff) || ( U) T^T^---TS%) = (u<8> e (/)>Ct«®e(s)>, 
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where we have again appealed to lemma [T] and theorem [3J Now note that 

\{u®4> k ,I$v®<p k ) - (u®e(f),R k v®e(g))\ 

<\(u®(ip k - e(/)), Rt v®ip k )\ + \(u® e(f), R k v® (<p k - e{g))) 

< \\u\\ \\v\\ (U k - e(/)|| H^ll + ||e(/)|| \\ Vk - e(g)\\ 0, 

where we have used that R k is unitary. Therefore 

(u ® e(f),R k v ® e(g)) k ^°°> (u ® e(/), C/ t u <g> e(g)) 

for any u, v £ Tt and simple functions f,g S L 2 ([0, TV]; C") with dyadic rational 
jump points. As the latter are dense in L 2 ([0, N]; C"), a similar approximation 
argument shows that R k converges weakly to Ut as k — > oo for every t £ [0, TV]. 

It remains to strengthen weak convergence to strong convergence uniformly on 
[0, N]. First, note that as R k and Ut are all unitary, weak convergence of R k to Ut 
(which is of course equivalent to weak convergence of (Rt )* to U£) already implies 
that (Rt)* — ► £/ t * strongly as k — ► oo for every fixed time t € [0,JV]. To prove that 
the convergence is in fact uniform, we will utilize another implication of theorem [3] 

It is convenient to extend the Fock space to two-sided time, i.e., we will consider 
the ampliations of all our operators to the extended Fock space T = T s (L 2 (K; C n ) ) = 
T-®J- ', where T- = T is the negative time portion of the two-sided Fock space. We 
now define the two-sided shift 6 t : L 2 (R; C") -> L 2 (R; C") as § t f(s) = f(t + s), and 
by Qt '■ J~ —* T its second quantization. Note that 0t is a strongly continuous one- 
parameter unitary group, and that the cocycle property reads Ut+ S = U s OlUtO s , 
etc., in terms of the two-sided shift. Now define on the two-sided Fock space the 
operators 

v t = 9 t u;, s k = e 2 - k (R k y. 

Then it is immediate from the cocycle property that Vt defines a strongly continu- 
ous unitary (hence contraction) semigroup on Tt ® T and that the unitary (hence 
contraction) S k is such that (S" fc ) L* 2 ^ J = Q 2 - k yt2 k \(R k )* ■ Moreover, for any dyadic 
rational t 

Wv-W 

for k sufficiently large, as 6 t is an isometry (here tp-®ip £ T- ®Tt® T = Tt® F). 
As vectors of the form ®tf) are total in Tt ®F and as we have already established 
that \\(R k )*ip — U£ip\\ — ► as k — > oo, we find that for every fixed dyadic rational t 

||(S fe )L t2fe J i/)-Vti/)\\ for £H®f. 

As Vt is strongly continuous and the dyadic rationals are dense, this evidently holds 
for every fixed t. It remains to apply the implication 2=>3 of theorem [3] □ 
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TheoremUl 2=>1. Fix a, (3 € C™, v G H, and t e [0,N], and let us choose the 
sequences ip k — e(a/[ 0! 2- fc ]) and <p — e {PI[a,2- k ])- Then we can estimate 



||(fl*tfV)L«%_i*/» w f = SU p \( u ,(RWv k )l t2h ) v - T^v)\ 2 

uETt 
II«II<1 

\{u®e{aI [0tN] ),(R k - U t )v®e{f3I [0 , N] ))\ 2 



||u||<l 

||(itf-fr t )tt<M/%AT])ll 2 

e \P\ 2 N 

2\\v® e(/3J [0;J v])|| 2 - 2Re((i? t fc « g eCftfajy]), t/« « ® e(pi [0iN] ))) 

e \l3\ 2 N 

But by assumption — U*) — > as k — > oo for all if> G Ji ® J 7 , so that in 

particular (R k v ® e(/3Jr 0l jvi), f/t u ® e {flI[o,N])) — * \\ v ® e (/3^[o,iV])l| 2 - We thus obtain 

IKiifc^V*)^*]^ _ -*=^ o for all « g ft, t > 0. 

The result now follows by appealing to the implication 2=>1 of theorem [3] □ 

4.2. Proof of theorem [2j The proof of theorem [5] is chiefly a matter of straight- 
forward computation. We make use of one simple trick: the trivial identities 

e ix = 1 + x f(x) = l + ix + x 2 g(x) 

and the conditions XjXo = 0, (xo, MjXo) — allow us to cancel those terms in the 
expression for 2 k {R k '^ v — I)u which diverge as k — > oo. 



Theorem^ Throughout the proof we fix a, ft G C™ and u G H., and we define 
ij) k = Ti k x k {a), ip k — ir k x k (P)- By theorem Q] it suffices to show that 

\\2 k (R k ^ k ^ k - I)u - ^ al3 u\\ 0, 

where S£ a ^u is given in lemma [2] in terms of the coefficients defined in section HOI 
For the time being, let us fix k. As H k is self-adjoint, we may assume by the 
spectral theorem that TC ® JC = L 2 (fl k ) for some measure space (f2 fc , S fc , P k ) and 
that H k acts on L 2 (fl k ) by pointwise multiplication (H k ip)(uj) = h k (ui)ip(ui) for all 
%jj G L 2 (£l k ). We represent the vectors u ® Xj ancl w <8> Xj m L 2 (Jl k ) as itj(w) and 
Vj(ui), respectively. As u®xo and v<8>Xo are in Dom(iJ fe ) (so that h k (lo)uo(lo) and 
ft fe (a;)t;o(co>) are square integrable), the trivial identity e lx = 1 + ix + x 2 g(x) gives 
that 

(v ® xo, 2 k {e lH " 2 ~ k -I)u® X a) = i ! v {u>)*h k {u)u Q {u) P k (duj) 



+ 2- k J {h k (cu)vo(u}))* g{h k {uj)2- k ) h k {uj)u {uj) P k (diu) 

= i{v®xo,H k u®xo)+2~ k {H k v®Xo,g{H k 2- k )H k u®Xo)- 
Using the trivial identity e lx = 1 + x f(x), we similarly obtain for p = 1, . . . , n 

(v ® Xp , 2 k (e mk2 ~ k -I)u® X o) = (v® x P , f(H k 2- k ) H k u ® X o), 
(v ® xo, 2 k (e lHk2 ~ k ~I)u®x P ) = (H k v ® xo, f(H k 2~ k ) u ® Xp ). 
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Using that XjXo = an d (xo, MjXo) = 0, a simple computation gives the following. 

(v® X fe («), 2 k {e^ k - I) u ® X fc (/3)> - 

r 

i^2(v <8 xo, Hju O Vjxo) 

3 = 1 

m 

+ (G j v® t j, j xo,9{H k 2- k )G f u®[i f xo) 

n m 

p=i j=i 

+ I]( G ^ ® MjXo, f(H k 2- k ) u ® x P ) /3p 
p=i j=i 

n 

+ ^ «; (« ® Xp, (e^ 2 " - /) « ® X 9 > /3p + 0(H 2- fe /2). 

p,g=l 

To proceed, let us define 

(n n n \ 

^ a£(JV M - + a * P M p + 51 L pPp + K ) u - 
p,g=l p=l p=l / 

Using the definition of N pq , M p , L p , K in section l2~4l we obtain 

r 

3=1 
m 

+ 51 {Gjv® ^jXo,g{F)Gj,u® fj,j>xo) 

n m 

+ 51 a p 53^ ^(^) g j u ^xo> 

p=l j=l 

n m 

p=i j=i 

n 

+ 5] o$ <« ® Xp, (e <p - /) u ® x 9 ) /V 

Thus evidently we obtain 

sup |(^x fc («),2 fc (e^ fc2 " -I) U ®x fe (/3)) ~ («,^tt)| -^0, 



provided that f[H k 2~ k ) -> /(F), g{H k 2~ k ) -> c/(F) and e^* 2 " -► e li? strongly 
as fc — » oo. This is indeed the case due to [2H1 theorems VIII. 20 and VIII. 25]. 
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To complete the proof, note that we can write 
|| 2 *(i2fc;*V -l) u -^u\\ = sup \{v,2 k (RW* k - I)u)-{v,^i 



r||<l 



sup 

.en 
IMKl 



2 k (v® X k (a),e* H ~ 2 ~ u® X k {P)) _ k , \ 

L \V,U) 



\x k (a)\\ \\x k (P)\\ 



< 



sn Pven,\\v\\<i \(v®X k {a),2 k (e 



i)u® x k m-{v,^ u)\ 



ix fc («)ii \\x k m\ 



k ( {X k {a),x k m 



1 11 + 



ix fe (")i! \\x k m\ ) iix fe («)ii \\x k m\ 

which converges to zero as k — > oo. The claim has been established. 
4.3. Proof of lemma [3j By our assumptions, we have 

<X fe («),xW fe = ((n k x k (a)r 2k ,(7r k x k (P)f 2k ) — 



□ 



for every a, ft £ C". For a, ft — we find that (xcb Xo) 2 '' - * 1 as k — » oo, so we must 
have |(X0) Xo)| = 1< But xo an d Xo are unit vectors, so xo = e^xo for some <j) £ K 
such that e*^ 2 — > 1 as k — ► oo. We may now compute for arbitrary a, ft £ C" 



<X»,xW=e^ 2 ' 



«ii=i 



— >oo 



exp 



<>3=1 



Substituting a h- ia and differentiating, wc find that 



dt 



t=o 



for every a,/? £ C". We thus find that (xj,Xj) = e ^ f° r every j, which implies 
that Xj = e Xj as Xj and Xj are unit vectors. The proof is complete. □ 

4.4. Proof of lemma [4l As the coefficients are bounded, it suffices to verify the 
Hudson-Parthasarathy conditions in remark [TJ We first show that 

n n 

NpjN q j — Spg, Nj p Nj q = S pq . 

3=1 3 = 1 

To see this, note that for all u,v £ Ti., the following identities hold true: 



(u,v)S pq = {u®x P ,e lF e lF v (g> x?) = (« ® X P , e lF Pe lF v® Xq ), 

-iFAF „,,<>,, \ _ /., o „, „-iF d„jF , 



(u, u) <S P9 = (u <g> x P , e ^ e li< u (g> x<?) = ( u ® Xp, e li< Pe l * u <g> x 9 ) 

for p, g = 1, . . . , n, where P := J2r=i T ® XrX* is the orthogonal projection onto 
Tt®S and we have used that XjS C 5 in the last step. Using the definition of N pq , 
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the claim is easily established. The next condition that we will check is 

n 

k + k* = -Y / l p l;. 

p =i 

Since 

n n n / n \ 

K + K* = G l (W l3 + W$G jt Yl L p L p = E G « E Y ? Y I* G *> 
the claim would follow immediately if we can show that for all i and j 

n 

p=l 

To see this, we begin by noting that —f(x)f(x)* = 2 (cos(x) — l)/x 2 . Furthermore, 
since 1-ijXo € S and XjS C S, we find that f(F)* v ® ^jXo &'H®S. Therefore 

f(F)f(F)* v ® MjX o - f(F)Pf(Fy v ® MjXo , 

and we can write 

/ 2cos(F) — 2/ \ - 

/utgi^xo, ^ v®l^jXo) = -(u® fiiXo, f(F)Pf(F)* v ® fijxo)- 

Using the definition of Yf, we find that for all w, u eW 

>T-/ \ / 2cos(F)-2J 

- 5j< U) 5f y/* v) = (u® tuxo, y 



p=i x 
On the other hand, note that g(x) + g{x)*~ = 2 (cos(a;) — l)/x 2 . We thus obtain 

, rx,*x x / 2 cost F)- 21 
(u, (Wij + W^) v) = ( u ® HiXo, p-2 v ® VjXo 

for all u, v G Ti from the definition of Wij. The claim is established. 
The last condition that we need to check reads 

n 

M p = -Y,N m L* q . 

9=1 

To see this, we begin by noting that —e tx f{x)* — f(x). Since fijXo & S and 
XjS C S, we find that f(F)* GjV ® fijXo <= <5 for all j = 1, . . . , to. Therefore 

/(F) ® wo = -e tp Pf(F)* GjV ® h jXo . 

We obtain 

E< M ® %» /(*0 G> ® HjXo) = - J2(u ® x P , e iF Pf(F)* G jV ® fi jXo ). 
i=i j=i 

Using the definitions of N pq , Y 9 , L q and M p , we find that the claim holds true. □ 
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